ON THE VIBRATIONS OF LUMPED PARAMETER SYSTEMS GOVERNED 
BY DIFFERENTIAL-ALGEBRAIC EQUATIONS 



S. DARBHA, K. B. NAKSHATRALA, AND K. R. RAJAGOPAL 

Abstract. In this paper, we consider the vibratory motions of lumped parameter systems wherein 
the components of the system cannot be described by constitutive expressions for the force in terms 
of appropriate kinematical quantities. Such physical systems reduce to a system of differential- 
algebraic equations, which invariably need to be solved numerically. To illustrate the issues with 
clarity, we consider a simple system in which the dashpot is assumed to contain a "Bingham" fluid 
for which one cannot describe the force in the dashpot as a function of the velocity. On the other 
hand, one can express the velocity as a function of the force. 
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1. INTRODUCTION 



q \ The traditional approach to obtaining the governing equations for the vibratory motions of a 

lumped parameter system of springs, masses and dashpots is to write down the balance of linear 
lO ■ 

j> ■ momentum for the system and to provide constitutive expressions for the forces in the springs 



and dashpots in terms of appropriate kinematical quantities. This leads to a differential equation 
for the motion which can be solved by a proper choice of the initial conditions for the system. 
However, in many instances, one is not in a position to provide the constitutive expressions for the 
forces in the components of the system in terms of the kinematical quantities, rather one finds that 
constitutive expressions can be provided for the appropriate kinematical quantity in terms of the 
forces. Two examples which immediately come to mind are a dashpot wherein the fluid behaves like 
a "Bingham fluid" or the dissipation due to Coulomb friction. In the case of a dashpot, wherein the 
fluid responds like a "Bingham fluid" (see Figure [1]), we immediately observe that the force is not a 
function of the velocity, but the velocity is a function of the force. (If by a fluid, one means a body 
that cannot resist shear stress, then the notion of a "Bingham fluid" is untenable as there cannot 
be a yield stress below which the fluid resists shear stress.) It is also possible that the response is 
such that the relationship between the velocity and the force is not invertiblc (sec Reference |13j). 
Similarly, it is possible that the energy storage mechanism is such that one cannot express the force 
in terms of the displacement but can only express the displacement in terms of the force. Such 
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a situation arises when one comes across an energy storage mechanism described by the lumped 



parameter system depicted in Figure 2(a) , which consists of a spring and an inextensible string in 
parallel. The relation between the force and displacement is that depicted in Figure |2(b)| This 



response can be expressed as the displacement being a function of the force as given in Figure |2(c) 



In general, the constituents of the lumped parameters may only be characterizable with implicit 
constitutive relations between the forces and kinematical quantities. Recently, Rajagopal [12] has 
discussed the general framework that arises when considering such lumped parameter systems. 
The problem reduces to the solution of a system of differential-algebraic equations. In this paper, 
we consider such a system of equations, for a rather simple lumped parameter system consisting 
of a mass, spring and dashpot, the fluid in the dashpot being characterized by the constitutive 
expressions for a "Bingham fluid" , with a view towards explaining the interesting features of such 
a differential-algebraic system of equations. This study is just an initial foray into an area which 
holds much promise in view of the lumped parameter system being much more complicated. For 
instance, as we mentioned earlier, the components of the lumped parameter system may be such 
that one can only provide implicit constitutive relations to describe them, with neither the force nor 
the kinematical quantity being expressible in terms of the other. It is possible that a constituent 
of the lumped parameter might have both an energy storage and an energy dissipation mechanism 
built into them, that is the lumped object cannot be further split into an energy storing (spring) 
and an energy dissipating (dashpot) mechanism. That is, the object that is being "lumped" might 
be viscoelastic. An example of such a situation presents itself when the spring-dashpot system 
reflects the representation of a Maxwell fluid [10] or one that corresponds to more complicated rate 
type of fluids. The system could then be comprised of various such lumped objects and more than 
one kinematical quantity and force might be given in terms of implicit equations (see |12j for a 
more detailed discussion of such lumped parameter systems). Also complicating matters further 
could be dissipative systems that are Coulomb-like. Such systems, even within the framework of 
simpler components are far from being well understood (see [151 [1] for a discussion concerning the 
mathematical difficulties associated with the study of such systems). 

In this short paper, we simplify the problem greatly in order to highlight the salient aspects of 
the problem under consideration. We study a simple mass, spring and dashpot system wherein the 
spring, to keep matters simple, is a linear spring and the dashpot is a "Bingham dashpot" whose 
response is given in Figure HI We obtain the differential-algebraic system corresponding to the 
motion of such a system and solve it numerically for different initial conditions. We find interesting 
results that are in keeping with physical expectations. For instance, in the case of free vibrations, 
it is possible that the system could come to rest in a position wherein the spring is stretched, i.e., 
the equilibrium solution is not the solution corresponding to the unstretched spring. Thus, the 




(a) (b) 

Figure 1. Bingham fluid model (a) force- velocity relationship (b) velocity as a 
function of the force 

equilibrium solution cannot be obtained by minimizing the "energy" associated with the system. 
This result is quite different from what one obtains in the case of a classical viscoelastic dashpot 
wherein the equilibrium state is the unstretched state of the spring. We also study the problem of 
the system being subject to different types of forced vibrations, the results being in keeping with 
physical expectations. 



2. PROBLEM STATEMENT AND GOVERNING EQUATIONS 

Consider a mass-spring-dashpot system as shown in Figure ([3]) . We shall denote the deflection of 
the spring from its unstretched position by x(t). We place an inertial frame of reference at the point 
occupied by the mass and hence, we may represent the displacement of the mass from this position 
by x(t). Let m denote the mass. Let Fd(t) and F s (t), respectively, denote the forces applied by the 
dashpot and the spring on the mass, and let F(t) denote the applied external force acting on the 
mass. An application of the balance of linear momentum yields the following equation: 

(1) x(t) = — (F(t) - F d (t) - F s (t)) , 

m 

where a superposed dot denotes the time derivative. In the first-order form, the above equation 
takes the following form 

(2a) v(t) = ~[F(t) - F d (t) - F s (t)}, 

m 

(2b) x(t) = v(t), 

where v(t) is the velocity of the mass. 

One must now specify a constitutive relationship for the spring force F s (t) and dashpot force 
Fd(t) in terms of the kinematical variables, usually x(t) and v(t), so that one can determine the 
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(c) Displacement as a function of force 
Figure 2. A spring with an inextensible string in parallel 
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Figure 3. A mass-spring-dashpot system 



response of the mass. Generally, these constitutive relationships usually take the form 



(3a) 
(3b) 



= a(x{t),F s {t)), 
= (3(v(t),F d (t)). 
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where a and /3 are appropriate functions. However, one could have a much more complicated 
implicit constitutive relationship than (|3a|) and (|3b|) . (For instance, the spring could be a viscoelas- 
tic solid though one invariably uses the term spring to denote a purely elastic body.) Then, the 
governing equations for the mass-spring-dashpot system may be written as 

(4a) v(t) = ±[F(t) - F.(t) - F d (t)], 

(4b) x{t) = v(t), 

(4c) = a(x{t),F s {t)), 

(4d) Q = p(v(t),F d (t)). 

Now the variables are x(t), v(t), F d {t) and F s (t). In addition to the above equations, we need to 
specify (consistent) initial conditions. In the next subsection, we shall consider specific models for 
the spring and dashpot. 

2.1. Specific models for spring and dashpot. In this paper we consider the case of spring, 
which is linear, and the dashpot is assumed to contain a classical Bingham fluid. The spring 
constant will be represented by k. The constitutive equation for the spring may be written as 

(5) x(t) = ±F 8 (t). 

To provide the constitutive relationship for the dashpot containing a Bingham fluid, we require two 
non-negative constants: /J, s Nf and 7, where [i s is coefficient of static friction, and Nf is the normal 
force on the mass. The constitutive relationship (I4dj) takes the following explicit form: 

(o) m = { mn<^"t 

[ ~l(Fd(t) - sgn[Ji(t)]»i.JV» \F d (t)\ > p.N, 
where 7 > is the slope and n s Nf > is the threshold force (see Figure EJ). For the case of 
> VsNf it is easy to see that 

(7a) sgn[v(*)] = sgn[F d (t)] 

(7b) 1^)1=7(1^(01-^/) 

Using equation © we eliminate x(t) in equations fljaj) - (j4d|) . and obtain the following equations: 

(8a) v(t) = ~[F(t) - F s (t) - F d (t)} 

m 

(8b) F s (t)=kv(t) 

(8c) = P(v(t),F d (t)) 

The above system of equations is a system of semi- explicit differential-algebraic equations (DAEs) 
(for example, see Reference [2]). The first two equations are ordinary differential equations (ODEs), 



x(t) = v(t) 
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Figure 4. A pictorial description of Bingham model 

and the third equation is an algebraic constraint. It is, in general, not possible to find analytical 
solutions for these system of equations, and one must resort to numerical solutions. In a subsequent 
section, we shall present a numerical algorithm to solve these differential-algebraic equations. Some 
representative references on differential-algebraic equations are [21 EJ El [TJ El [9] . 

Remark 2.1. When 7 — > 00, Figure^ resembles the response for Coulomb damping. However, 
it should be borne in mind that one cannot obtain the frictional response of solids by a limiting 
process for the frictional response of fluids. While the response curves might look the same, the 
philosophical underpinnings as well as other relevant physical issues are quite different. 

2.2. Computation of a solution for the differential-algebraic system of equations. Clas- 
sical textbooks on vibrations such as the one written by Meirovitch supposedly provide a method 
for obtaining the response of the spring-mass system subject to Coulomb damping (see Meirovitch 
[111 pages 31-34]). The method consists of piecing together solutions corresponding to the regimes 
when x > and x < 0. The resulting solution obtained by Meirovitch is smooth, as can be seen 
from the figure on page 34 of |llj . The procedure adopted by Meirovitch runs into trouble when 
one considers the following initial conditions: x(0) = xq, where x(0) = and \xq\ < f . The 
main difficulty lies in the inability to specify when x(0) = 0. 

Put differently, the main error of Meirovitch's analysis stems from the fact that when the initial 
conditions are such that |xo| < ^ s Nf/k and x(0) = 0, one cannot use the differential equation (JTJ) to 
describe the system as F^ is undefined, within the framework that Meirovitch used to describe the 
problem. The differential equation ([I]) is incapable of describing the response of the system. One 
can proceed in either one of two ways. Use the differential-algebraic system (P), ([5]) and ([6]). This 
framework does allow one to describe the system as we now do not choose to provide a constitutive 
relation for Fj, which we are not in a position to provide, rather we choose to treat the force F^ 
as an unknown that is obtained as a solution to the system (pTJ) , (0) and (0). The other way of 
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studying the problem that is well-posed is to describe the physical problem within the context of 
a differential inclusion, as considered by Filippov [5]. The approach adopted by Filippov [5] is 
described briefly in the discussion that follows. We shall, however, consider the problem within 
the context of a differential-algebraic system in some detail as such an approach can be adopted 
to study a whole host of problems wherein the forces in the springs and the dashpots are given in 
terms of implicit constitutive relations that relate the forces to appropriate kinematical variables. 

If one is to interpret the frictional force in terms of the signum function, with sgn[0] = 0, then 
one does not have a classical solution to the problem. One may appeal to intuition in this case and 
claim that the solution is x(t) = xq and x(t) = 0. However, this solution, as we will see, is not a 
solution in the classical sense. Neither does Meirovitch specify the sense in which the function he 
has constructed actually solves the problem. The absence of the specification of the sense to which 
one has a solution to the problem, namely the class of functions to which the solution belongs, can 
lend itself to misinterpretation of the function being a solution in the classical sense. Meirovitch 
takes an intuitive engineering approach to the problem and does not discuss the mathematical 
subtleties of the solution. 

2.2.1. Non-existence of classical solutions. Before we proceed further, it is better to define what we 
mean by a classical solution or a solution in the sense of Caratheodory [4j. Let z(t) € Si n represent 
the state of a dynamical system at time t and f(t,z(t)) denote a function mapping 3? x K n —> W 1 . 
Let zq represent the state of the dynamical system at t = 0. 

Definition 2.2. Consider an n th degree ordinary differential equation with the following represen- 
tation: 

(9) z(t) = f(t,z(t)), z(0) = z . 

We will say that <p(t) is a solution in the sense of Caratheodory to the differential equation if 

• (f)(0) = Zq, 

• (f)(t) is absolutely continuous, and 

• t(t)=Zo + J*f(T,<l>(T))dT. 

We will define <f>(t) to be an equilibrium solution if, in addition to it being a solution, (f)(t) = (f)(0) 
for all t > 0. 

The non-existence of classical solutions stems from the nature of dry friction. Physically, there is 
a continuum of equilibria which correspond to the velocity of the mass being equal to zero and the 
deflection of the spring not being sufficiently large for the spring force to overcome the friction force 
and cause motion. At all these equilibria, the acceleration and velocity of the mass is zero. However, 

if we were to set the acceleration to be zero when the velocity is zero, we see that the friction force 
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must balance the spring force and hence, the friction force cannot be specified independently of 
the initial deflection of the spring. In fact, depending on the initial deflection of the spring, the 
Coulomb damping force can take any value in (— fj, s Nf , fJ> s Nf) when the velocity is zero. Hence, the 
specification of Coulomb damping force, Fj, cannot be performed solely on the basis of velocity as 
it is dependent on the deflection of the spring when x = 0. 

If one were to describe the damping force to be = —fj, s Nf sgn[x], the implication is that the 
damping force is specified when the velocity is zero and this is not physically the case. It is exactly 
this difficulty that leads to the non-existence of classical solutions for almost all initial conditions, 



where the initial velocity is zero and the initial deflection of the spring is no more than — ^ 
magnitude. 



m 



2.2.2. Existence of solutions in the sense of Filippov. Let us begin with the definition of a solution 
in the sense of Filippov |14j . We will use the following notation: For every z € 3? n , the term B(z, S) 
represents the set of all points in ffl 1 which are no farther than 5 from z. For every set N C 3i n , let 
fj.(N) denote the measure of the set. Let A, B be two subsets of 3? n ; their difference, represented 
by A — B indicates the set of points in A that do not belong to B. 

In Reference [5], Filippov considered the following differential inclusion: 



(10) 



x £ F(x,t) 



where F(x, t) is the possible set of values that x can take. Filippov indicated that if for every (x, t), 
F(x,t) is non-empty and convex, and the distance between the sets F(x ,t ) and F(x,t) tends to 
zero as x — > x and t — > t , then a solution (that is absolutely continuous and satisfies the differential 
inclusion (jlOp almost everywhere in time) exists. 

Consider the set of values the right-hand side of equations ()8ap and (|8b|) can take, which can be 
written as 



(11a) 
(lib) 
(11c) 



F + 



F- 



F n 



F(t) 



F(t) 



kv 



kv 



for v(t) > 
for v(t) < 



X[ F (t)-F a -F d ] 




, for v(t) = 0, where \F^\ < n s Nf 



For the case v ^ 0, the set F(x, t) is a singleton set. For the case v = 0, F(x, t) is a convex set and 
the other conditions are met. Solutions exist and are guaranteed to pick appropriate value of F$ 
from the set {\F^\ < fj, s Nf} when the velocity is zero. 



2.3. A numerical algorithm for solving a system of differential-algebraic equations. 

It is well-known that Backward Difference Formulae (BDF) are stable and accurate (and hence 
more suitable) for solving DAEs numerically [2j [3] . The backward Euler scheme is the simplest 
member of BDF. We shall now discretize equations ([5]) using the backward Euler time stepping 
scheme. To this end, we shall discretize the time interval of interest into A?" time instants denoted 
by t n {n = 0, • • ■ ,N). For simplicity, we shall assume uniform time steps, and shall denote it by 
At := t n — t n -\. We shall denote the (time) discretized version of a given quantity z(t) at the 
instant of time t n as 

(12) z {n) « z(t = t n ) n = 0,---,N 

The corresponding discretized equations using the backward Euler scheme at the instant of time 
t n +l can De written as 

(13a) ± k> +1 ) - «W] = 1 [F^ - Fj n+1 ) - i> +1) ] 

(13b) Kt ^ n+1) ~ F ° n ^ = k v(n+1) 

(13c) = /3(*/ n+1 \i> +1) ) 

Using equation ()13b|) we can write i ? i n+1 ^ in terms of the variable v( n+1 ) as 
(14) Fj n+1 ) = + At k v in+1) 

Using the above equation, equation (|13a|) can be written as 



(15a) 



1 + At 2 — 

m 



y (n+l) = p{n+l) _ p (n+l) 

m V d 



(15b) := + - F,W 

From equation (|15aj) (and noting equation ([7a|) for the case of > fJ-sNf) we have the following 

useful result 

(16) sgn[F( n+1 )] = sgn[Fj n+1) ] 

which is valid for all cases (that is, for both |Fj n+1 ^| < [i s Nf and l-F^i™ | > fJ> s Nf). The above 
result can be utilized quite effectively in a numerical scheme as F( n+1 ) is known whereas F^ n+1 ^ is 
not known as it is part of the solution for the instant of time t n +l- We shall employ a predictor- 
corrector-type scheme to calculate the values for the instant of time t n +i assuming that all the 
values at the prior instants of time are known. The scheme is outlined in Algorithm [TJ 

Remark 2.3. In the case of the standard Bingham model, the (discrete) solution from t n to t n+ \ 
can be solved analytically. But for more complicated models, one may have to employ a Newton- 
Raphson scheme to solve the nonlinear equation, which can be written just in terms of i^j n+1 ) . 



Algorithm 1 Predictor-corrector algorithm for advancing the solution from t n to i n +i 

1: Input: k, m, At, fj, s N f , v^ n \ F s (n) , F^ n+1 ^ 
2: Output: x( n+1 \ v^, F { d n+1 \ F s (n+1) 

3: Predictor step: Calculate the predictor := + F( n+1 ) - F^ n) 

4: Corrector step: 

5: if |F( n+1 )| < n s N f then 

6: = and Fj n+1) = 

7: else 



,(n+l) 



t 2i| 7 At 



7[1+At : 

,,+i) = 7 (Fj" +1) - sgn[Fi" +1 VsiV/) 



7 (1 + At 2 A) sgn^^+^j^iV/ + 



9: U 
10: end if 

11: Calculate Fj n+1 ' ) and x( n+1 ) using 



F jn + 1) = pin) +Afk v (n+l) ^ x (n+l) = * F (n+l) 



3. NUMERICAL EXAMPLES 
In all the numerical examples we shall take 7 = 1, /^ s iV/ = 1, m = 1, k = 100, and At = 10 -4 . 

3.1. Non-zero external forcing function. We shall take x(0) = and Fd(0) = 0. These 
conditions will imply that the corresponding (consistent) initial conditions for the other variables 
will be v(0) =0 (from equation ([6])) and F s (0) = (from equation ©). Two different external 
forcing functions are considered 

0.5sm(5vrf) < t < 1, 
t > 1, 

f 10.0sin(5vrt) 0, < i < 1 

17b F 2 t = 

[ t > 1. 

The numerical results for the forcing function F\(t) are presented in Figure The numerical results 
for the forcing function F2(t) are presented in Figures [6] - [9l 

Remark 3.1. The forcing functions are chosen in such a way that 

(18) max[Fi(t)] < fi s N f and max[F 2 (t)] > [i 8 Nf 
The analytical solution for the first case is 

(19) F d (t) = Fi(t), F 8 (t) = 0, x(t) = 0, v(t) = 
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(17a) 




Figure 5. Non-zero forcing function using equation (|17ap . 




Figure 6. Displacement x(t) for non-zero forcing function given in equation (|17b|) . 



As one can see from Figure \T7a\ the numerical solution matches the analytical solution quite well. 



3.2. Non-zero initial displacement. We shall take the external force F(t) = 0, and take Fd(0) = 

0, which implies that the consistent initial condition for the velocity will be v(0) =0 (from equation 

([6])). Two different initial displacements are considered: x(t = 0) = 0.005 (which implies that 

F s (t = 0) = 0.5), and x(t = 0) = 0.5 (which implies that F s (t = 0) = 50). The numerical results 

li 



1 CI v l l l l l I 

0.5 1 1.5 2 2.5 3 

time 



Figure 7. Velocity v(t) using non-zero forcing function given in equation (j!7b[) . 




0.5 1 1.5 2 2.5 3 

time 

Figure 8. Force in the spring F s (t) using non-zero forcing function given in equa- 
tion (|17b|) . 

corresponding to x(0) = 0.005 are presented in Figures [TU1 The numerical results for the case 
x(0) = 0.5 are presented in Figures [TT1 - [T4l 

4. CONCLUSIONS 

In this short study we have considered problems of vibratory motion of lumped parameter systems 
wherein the forces in the springs and dashpots belonging to the system cannot be given as functions 
of the kinematical quantities; rather the kinematical quantities are expressed as functions of the 
forces. This leads to the system being described by a set of differential-algebraic equations. We 
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Figure 9. Force in the dashpot Fd(t) using non-zero forcing function given in equa- 
tion (|17bj) . 
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Figure 10. Non-zero initial displacement of x(0) = 0.005. 

illustrate the type of equations one has to deal with by considering a simple system, for the sake of 
simplicity of illustration. We considered a very simple mass-spring-dashpot configuration with the 
spring being a linear spring and a dashpot that contains a Bingham fluid. We subjected this system 
to a class of forcing functions and initial conditions and found the solutions to be in keeping with 
physical expectations. In the future, we plan to study the vibration of systems whose components 
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Figure 11. Displacement x(t) using non-zero initial displacement of x(0) = 0.5. 




Figure 12. Velocity v(t) using non-zero initial displacement of x(0) = 0.5. 

are described by implicit constitutive relations between the forces and the kinematical variables 
which cannot be simplified either to the forces being given in terms of the kinematics or vice- versa. 
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